Finite band inversion of ARPES in Bi 2 Sr 2 CaCu 2 08+5 and comparison with optics 
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Using a maximum entropy technique within a finite band Eliashberg formalism we extract from 
recent high accuracy nodal direction angular resolved photo-emission spectroscopy (ARPES) data in 
optimally doped Bi2Sr2CaCu208+i (Bi2212) a quasiparticle electron-boson spectral density. Both 
normal and superconducting state with d-wave gap symmetry are treated. Finite and infinite band 
results are considered and contrasted. We compare with results obtained for the related transport 
spectral density which follows from a similar inversion of optical data. We discuss the implication 
of our results for quasiparticle renormalizations in the antinodal direction. 

PACS numbers: 74.20.Mn 74.25.Gz 74.72.-h 



I. INTRODUCTION 

Maximum entropy techniques have proved useful in at- 
tempts to extract boson information either from angular 
resolved photo emission spectroscopy (ARPES) data or 
from the optical conductivityi^^^ What is recovered 
is the electron-boson spectral density, 7 2 x(w), which de- 
scribes the effective interaction between two electrons due 
to the exchange of a boson which could be a phonon, a 
spin fluctuation, or some other excitation such as a plas- 
mon. In this way one can obtain some insight into the 
nature of this interaction and by implication about mech- 
anism of superconductivity. In particular, coupling to a 
resonance peak in 7 2 x(cj) can lead to peaks or kinks in 
measured quantities. As we will see, even a structure- 
less background can be picked up in our inversion pro- 
cess. Other approaches to the analysis of such structures 
have also been applied, see for example Mishchenko and 
Nagaosa 6 who employed the t — J-model. 

ARPES gives information on the quasiparticle self en- 
ergy while optical data can be expressed in terms of an 
optical self energy which is related to, but is different 
from the quasiparticle self energy. There are two main 
differences. The first is that optics involves the current- 
current correlation function which can be expressed in 
terms of a two particle Green's function and there can be 
vertex corrections while ARPES requires only the knowl- 
edge of the one-particle Green's function. Optics deals 
with a transport process and transport lifetimes are not 
the same as quasiparticle lifetimes. For example, the 
effectiveness of a scattering process in depleting a cur- 
rent depends strongly on the final state, i.e.: backward 
as opposed to forward scattering. Secondly, optics deals 
with a momentum average while ARPES is momentum 
specific. There are other complications: ARPES mea- 
sures renormalized quasiparticle energies directly and to 
extract from this data the quasiparticle self energy it is 
necessary to know the bare dispersion relation. Often, for 
energies not too far from the chemical potential, a linear 



dispersion is assumed and its slope is determined from 
an assumption that the self energy crosses zero at some 
energy ~ 400 meV in the recent high accuracy data in 
Bi 2 Sr2CaCu208+,5 (Bi22I2). 7 A similar ambiguity arises 
in optics: to get the optical self energy from the reflec- 
tivity data it is necessary to specify a value for the di- 
electric constant e at infinity. Also the plasma frequency 
is needed and this quantity is not so well known in the 
cuprates. In some determination it requires an assump- 
tion about where the band of interest ends as there are 
overlaps with higher bands, creating ambiguity. 

Certainly we are dealing with a finite band situation. 
For a simple first neighbor only tight binding band with 
hoping t, the band width W — 8t. Estimates based on 
tight binding fits to local density approximation (LDA) 
band structure calculations give values of the order of 
350 — 450 meV, Ref. Q for t while fits to experiment can 
give somewhat smaller values of order 200 meV. In all 
cases, of course, higher nearest neighbor hoping is also 
present. 

So far maximum entropy inversion techniques have in- 
volved using as an effective low energy theory for the self 
energy, the Eliashberg equation generalized to include 
any boson-exchange mechanism and not just phonons. 
But these are written for infinite bands. A deficiency of 
such an approach is that the self energy cannot change 
sign in an infinite band nor can its optical counterpart. 
On the other hand in finite band formulations a change 
of sign occurs naturally and is a robust feature of the 
formulationi2ii£ In early ARPES experiments the real 
part of the quasiparticle self energy was simply assumed 
to go to zero at ~ 300 meV. In the work of Meevasana 
et a/Ai for the La2- a; Sr 2 ;Cu04 (LSCO) series the renor- 
malized dispersions were found to cross the bare LDA 
bands at most dopings considered, with the energy of the 
crossing falling roughly in the range of 400 to 600 meV. 
Further, in the Bi2212 series of Hwang et al^ the op- 
tical self energy was also found to go through zero. On 
the theoretical side Cappelutti and PietroneroAS, noted 
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that in an electron-phonon model the real part of the self 
energy always goes through a zero at some finite energy. 
Later Knigavko and Carbotte^ established that for cou- 
pling to an Einstein mode the zero occurs at ~ ^JlueW/2 
where lue is the frequency of the oscillator. For optics the 
crossing is at higher energies, ~ \[2 times its quasiparti- 
cle counterpart to logarithmic accuracy. 

In this paper we start by generalizing the maximum 
entropy inversion technique for the quasiparticle self en- 
ergy used by Shi et al£ to finite bands. We use this 
new formalism to study how the value of the band width 
W changes results obtained for the electron-boson spec- 
tral density I 2 x(cu)- First we consider the normal state 
and generalize the procedure later on to deal with su- 
perconductivity. For this purpose it is necessary to have 
finite W Eliashberg equations with d-wave symmetry for 
the superconducting gap. We use these equations to de- 
rive from the data of Zhang et alZ the spectral functions 
I 2 x(ou) at low temperatures in the superconducting state. 
The available data is for the nodal direction only. In prin- 
ciple, the spectral density depends on direction and so 
antinodal results would be expected to be different. On 
the other hand, optics involves a transport spectral den- 
sity which is an average over all directions in momentum 
space. While this means that ARPES and optics cannot 
be compared directly, 13 wc, nevertheless, use the spec- 
tral density I 2 x{cu)tr taken from the optical literature to 
get some approximate information on a quantity which 
should be close to the angular averaged quasiparticle self 
energy. We comment on the points of agreement as well 
as the disagreements that are found. 

In Sec. |TT] we provide details of our formalism and de- 
scribe results for the normal state. Section Unl deals with 
the superconducting state and also provides a compari- 
son with optics. Finally, Sec. IIVI gives a brief summary 
and conclusions. Mathematical details are found in Ap- 
pendix [XJ 



Eliashberg theory in the normal state- 1 
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Here, N(u>') is the fully renormalized electronic density of 
states and carries the information on finite band effects, 
f(oj) and n(v) are the Fermi and Bose distribution func- 
tions, and d is an infinitesimal positive parameter. For 
finite bands N(cu) = Nq(Q) and K(lu + id, v) reduces to a 
closed form. Here No(G) is the bare electronic density of 
states at the Fermi energy taken to be a constant. 

It is clear that both real and imaginary part of E(w-H<5) 
are related to the desired spectral function through a con- 
volution integral to which maximum entropy techniques 
apply and either can be used. Inversion of self energy 
data on H{oj+id) in finite bands requires additional infor- 
mation on N(lu). This complication, however, can easily 
be handled. The quasiparticle spectral density A(k, lu) is 
related to the one-particle Green's function G(k, lu) by 



A(k, lu + id) = 3m {G(k, lu + id)} , 

7T 



with Dyson's equation 
G(k,ou + i5) 
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lu + id — £k — ^(lu + id) ' 
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where is the bare electron dispersion relation. For 
isotropic bands we can use e = e(k) to label the states 
instead of momentum and take the simplest finite band 
model for the bare density of states N (e) — 1/W for 
e in the interval [— W/2,W/2] and zero otherwise. The 
renormalized quasiparticle density of states is then 



N{lu) = / deN (e)A(e,Lu) 



II. FORMALISM, NORMAL STATE 

Maximum entropy techniques can be used to extract 
a spectral function I 2 x(cu) from the knowledge of the 
quasiparticle self energy 4- iS) related in integral 
form through a known kernel K{lu, v) (specified below), 
namely 



E(w + iS) 



dv K(lu + iS, v)I 2 x(v), 



(1) 



where / 2 x(w) is the electron-boson spectral density which 
describes the interaction of two electrons by the exchange 
of a boson of energy v. The kernel K(lu, v) is within 
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Let us assume we know T,(lu + iS) by some means, then 
N(lu) is known from Eq. © and the kernel Eq. Q for 
the maximum entropy inversion is now definite and the 
procedure can be carried out. 

ARPES experiments usually provide information only 
on the real part of the quasiparticle self energy, T,i(lu). 
Thus, we separate Eq. |T]) into its real and imaginary part 
and apply the maximum entropy method to deconvolutc 
only 



Si(w) 



dvK{Lu,v)I 2 X {u), 



(0) 
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with the kernel 



K{w,v) 
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(7) 

Here, V indicates that a principle part integral is to be 
taken. Eq. ([7]) differs from the one suggested by Shi et 
al£ in their inversion work in two regards. It contains the 
fully renormalized density of states N(u>) which accounts 
for finite band effects and the Bose distribution function 
n(v) is included because we do not want to be restricted 
to the low temperature range. 

If the imaginary part of the self energy, £2(0;), is not 
known from experimental data Eq. ([6]) cannot be applied 
directly to calculate N(u) and it is, therefore, required to 
develop an iterative, self consistent, formalism which will 
in the end allow us to extract the desired spectral func- 
tion J 2 x(w) in the finite band case from £i(w) alone. As 
the kernel ([7]) is based on Eliashberg theory, it is only 
natural to use the set of finite band d-wave Eliashberg 
equations ; 20 ! 21 ' 22 ! 23 as they are given in Appendix [A] to 
calculate £(w) for a given temperature T and a given 
spectral function 7 2 x(w). When this is done, the nor- 
malized quasiparticle density of states N(ui) of Eq. ^ 
can be evaluated for any choice of spectral density. The 
following self consistent procedure can be established: (1) 
An assumption is made for N(lo) and the simplest one, 
namely N(lo)/N (0) = 1/W in the interval [-W/2, W/2] 
will suffice. Here, W is the band width. (2) Eq. ([6]) is de- 
convolved using maximum entropy techniques and the 
experimental data on £i(w) for a given temperature T. 
In this step it is necessary to adjust W, which is an exter- 
nal parameter to the deconvolution, as is N(u>), for best 
data reproduction. The result is a first approximation to 
the desired spectral function I 2 x(u). (3) A solution of 
the finite band e?-wave Eliashberg equations (IA1|) based 
on this approximate function I 2 x(oj), the assumed value 
of W and the given temperature T is generated. From 
this solution the complex quasiparticle self energy £(w) is 
easily calculated and Eq. ([6]) can be solved to give a new 
fully renormalized density of states N(oj) and the proce- 
dure returns to step (1) until self consistency is reached. 
Another possibility is to start with step (1) and (2) from 
above. The approximate solution for I 2 x(u) can then 
be parameterized and a least squares fit procedure using 
the finite band d-wave Eliashberg equations (|A1[) can be 
employed to get the best fit to the experimental £i(w) 
even when the data is taken in the superconducting state. 
This is our preferred method. 

To get some understanding of how important finite 
band effects might be in maximum entropy inversions of 
the quasiparticle self energy data we proceed as follows. 
Hwang et al- have obtained from optical conductivity 
data results for the average transport spectral function 
7 2 x(w)tr in Bi2212 as a function of temperature and dop- 
ing. The solid curve of Fig. [1] reproduces their results 
for an overdoped sample (T c = 82 K, labeled BI82B) at 
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FIG. 1: The electron-boson spectral density I 2 \(lo) recovered 
from inversion of the real part of the quasiparticle self energy 
calculated with common I 2 x(w) (solid curve) but with differ- 
ent band width W = 00 (dashed), W = 2.5 eV (dotted), and 
W — 1.25 eV (dash-dotted). In all cases the inversion was 
carried out assuming an infinite band. 



T = 26 K. This function is used in Eqs. © and (J7J with 
W — > 00 so that N(u')/N (Q) = 1, to obtain the real 
part of the quasiparticle self energy of Eq. (JTJ) . [We will 
refer to it as the input / 2 x(o>).] Next, the procedure is 
reversed and maximum entropy techniques are used to 
derive from this numerical data a new spectral function 
J 2 x(w) which is shown as the dashed curve in Fig. [1] and 
which is seen to be almost identical to the input curve as 
it must be. Next we compute again £1(0;) from the same 
input 7 2 x(w) but now we apply a finite band cut off to 
the normal state version of the Eliashberg equations ()A1[) 
to generate a new set of numerical £i(w) data. This new 
set is again used as input for a maximum entropy infi- 
nite band deconvolution of Eq. ([6]) to yield a new model 
J 2 x(w). Results of this procedure for W — 2.5 eV (dotted 
curve) and W — 1.25 eV (dash-dotted curve) are given in 
Fig. Q] It is seen that the application of a finite band 
cutoff to the calculation of £i(w) has a major affect at 
energies beyond 100 meV where the new I 2 x(co) is con- 
siderably reduced over its input value. This was to be ex- 
pected since inverting in an infinite band model does not 
account for the reduction in self energy that is brought 
about by the decay in the effective electronic density of 
states around the bare band edge and beyond. It is clear 
that to do realistic inversions in the cuprates finite band 
effects need be accounted for. This is also required if the 
real part of the quasiparticle self energy is to cross zero 
at some finite energy, say 400 meV as was assumed in 
the analysis of Zhang et alZ experimental data. On the 
other hand, if one is mainly interested in the small energy 
region much less than the zero crossing at ~ 400 meV, fi- 
nite band effects make little difference except for lowering 
somewhat the peak around 60 meV in comparison to our 
input spectral function. 
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FIG. 2: Result for the electron-boson spectral density I 2 x(uj) 
obtained from inversion of Zhang et al. 7 - ARPES nodal di- 
rection Bi2212 data at T = 99 K. The dotted curve is for an 
infinite band, the green one for a band width W = 2eV, and 
the solid one for W — 1.2 eV. Note the reduction of spectral 
weight beyond ~ 100 meV as W is increased. Values of the 
mass enhancement factor are 0.81, 0.752, and 0.714, respec- 
tively. 



Having established our inversion technique, we next go 
to experiments. We started with the T — 99 K data from 
Fig. 4a of Ref. [7| and inverted it according to Eqs. © 
and ([7]) to recover the electron-boson spectral function 
I 2 x(w) in the specific case of the nodal direction at 99 K 
in the normal state. Results are presented in Fig. [2] Two 
different values of the band width, namely W = 1.2 eV 
(solid curve) and 2.0 eV (dashed curve) where used to- 
gether with W = oo (dotted curve). In all cases the 
same experimental data appears on the left hand side of 
Eq. © . For the dashed and dotted curves a band width 
W = 2 eV and W — oo, respectively, was imposed on the 
inversion procedure from outside. For the solid curve, 
the parameter W was allowed to vary. Instead, a con- 
straint, namely that the resulting self energy Xi(a->) be 
zero exactly at w = 400 meV was applied. This resulted 
in a value W = 1.2 eV. Of course, it needs to be rec- 
ognized that the ARPES experiments themselves do not 
tell us where the zero in £i(w) occurs. Some assumption 
on the bare dispersion is needed and the value 400 meV 
while respected in our inversions for the solid curve is, 
therefore, model dependent. If one had an independent, 
reliable estimate of the band width then this value could 
be used as a constraint in the inversion and this would 
yield an estimate of the energy at which renormalized 
and bare band dispersions meet. Returning to Fig. [2] it 
should now be clear why for a fixed set of Si(w) data in- 
creasing the value of W leads to smaller values of I 2 x{uj) 
at higher energies. Renormalization effects in this re- 
gion can be reduced due to a smaller value of I 2 x(u}) 
which, in turn, results in a reduced quasiparticle density 
of states. Finally, we note that the mass renormalization 



i \ Inversion from Ref. [6], T = 99K 

i i Inversion from Ref. [5], T = 102K 
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FIG. 3: Top frame: Electron-boson spectral density I 2 x(&) 
from finite band inversion (W = 1.2 eV) of ARPES data in 
nodal direction (Ref. [3]) in nodal direction in the normal state 
at T = 99 K (solid curve) . The dashed curve is for comparison 
and was obtained previously by the inversion of optical data 
(Ref. [H) and scaled down by a factor of 0.44 so as to account 
for differences between quasiparticle and transport quantities. 
Bottom frame: Real part of the quasiparticle energy at 
T — 99 K in the normal state. The open circles are the data 
of Zhang et al£ as read off their Fig. 4a. The solid curve rep- 
resents the result of the maximum entropy inversion of the 
data assuming a band width of W — 1.2 eV so as to get a 
zero in the self energy at ~ 400 meV as in the experimental 
data. The dashed curve was obtained using the I 2 x(&) shown 
by a dashed line in the top frame. A band with of 1.6 eV was 
chosen to, again, give a zero in £r(w) at ~ 400 meV. 



value A = 2 dui I 2 x{u)/uj is of the order 0.7 to 0.8 for 
optimally doped Bi2212. These values are considerably 
smaller than those determined from optics as we will dis- 
cuss later. There appears to be a factor of two difference 
in the magnitude between quasiparticle and transport 
electron-boson spectral density. 

Our best fit for I 2 \{uj) to the data of Ref. 0] is repro- 
duced as the solid curve in the top frame of Fig.[3]where it 
is compared with data on I 2 x(uj)t-r (dashed curve) which 
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was obtained from maximum entropy inversions of opti- 
cal data by Hwang et al.^ for a similar optimally doped 
Bi2212 sample. A scaling factor of 0.44 was applied to 
J 2 %(a;)tr in this case. The resulting mass enhancement 
A = 1.09 which is to be compared with the value of 0.72 
we found from ARPES. The above scaling factor was 
determined in an attempt to get the best possible agree- 
ment to the ARPES quasiparticle self energy Ei(u;) with- 
out changing the shape of / 2 x(w) tr . The corresponding 
self energy is shown as a dashed line in the bottom frame 
of Fig. [3] The solid curve derived from ARPES data, 
of course, fits data very well within the considered en- 
ergy range of [— 0.4eV, 0] for the experimental data of 
Ref. which are indicated by open circles. The dashed 
curve calculated from the rescaled I 2 x(^)tr (dashed line 
in the top frame of Fig. [3]) shows remarkable similar- 
ity. The main difference is due to the fact that the peak 
in 7 2 x(w)tr from optics at 62meV is stronger than the 
one in the ARPES data which is at 68 meV. This is not 
unexpected since optics produces an electron-boson spec- 
tral density which is an average over all momenta while 
ARPES is momentum specific, namely k is in the nodal 
direction. If we associate the peak with the interaction 
of the charge carriers with some spin fluctuations peak 
around (tt, tt) in the two-dimensional CuO Brillouin zone 
then we would expect this peak to be larger for scatter- 
ing in the antinodal direction and, therefore, larger in the 
average function of optics than is the nodal function of 
ARPES. 

The good agreement found here between the ARPES 
and optics derived spectral density noted in the top frame 
of Fig. [3] shows that both methods agree that there is 
strong coupling to an excitation at u> = 60 meV as well 
as a high energy background which extends to 400 meV. 
This cannot be due to phonons but finds a natural inter- 
pretation as coupling to spin fluctuations. 



III. SUPERCONDUCTING STATE AND 
RESULTS BASED ON OPTICS 

In the superconducting state the basic inversion pro- 
cedure outlined in the pervious section is no longer ap- 
plicable because Eqs. ^B) and ([2]) no longer apply. We, 
therefore, adopted the least squares fit procedure out- 
lined in connection with Eqs. ^ and (0 to perform the 
inversion of the nodal direction Ei(w) ARPES data of 
Zhang et alX In Fig. 3] we present our results of max- 
imum entropy inversion for temperatures T = 128 and 
99 K in the normal state and 35 and 17 K in the super- 
conducting state. In all cases we get a first estimate for 
I 2 x(w) from a deconvolution of Eq. ([6]) using the maxi- 
mum entropy method. This initial form is parameterized 
and then a least squares fit like procedure is applied to 
fit the theoretical values found from finite band 

Eliashberg theory, Eqs. (|A1|) . to experiment. The band 
width is kept at W = 1.2 eV. All curves show a peak 
around ~ 65 meV which is most prominent at 17 K in the 
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FIG. 4: Results for the electron-boson spectral density I 2 x(u) 
obtained by inversion of ARPES data by Zhang et ali^ along 
the nodal direction at different temperatures, namely T — 
128 K (dash-dotted line, A = 0.73), T = 99 K (dotted line, 
A = 0.72), T = 45 K (dashed line, A = 0.93, superconducting 
state), and T = 17 K (solid line, A = 1.12, superconducting 
state). A finite band width W = 1.2 eV was applied. 

superconducting state. As T increases this peak broad- 
ens somewhat and shifts towards higher energies. Besides 
this resonance like peak there is a large, structured back- 
ground which exists up to 400 meV. It consists of a valley 
with its lowest point around ~ 115 meV and additional 
structure beginning at energies ~ 150 meV. The real part 
of the self energy Si (uj) obtained from these spectra after 
solution of the Eliashberg equations (|A1|) is shown in the 
top frame of Fig. For clarity we do not show data but 
in all cases a tight fit was obtained, comparable in qual- 
ity to the fit shown in the bottom frame of Fig. [3] As it 
is of considerable interest to compare these results with 
optics we present the results of additional calculations in 
the bottom frame of Fig. [51 Here we used, as in the bot- 
tom frame of Fig. [31 for / 2 x(w) the spectra derived from 
the optical data reported by Hwang et al£ for a number 
of temperatures, namely 300 K, 200 K, 102 K, all in the 
normal state, and 72 K and 27 K in the superconduct- 
ing state. [The normal state results for T = 27 K (solid 
gray line) have also been included for comparison. In all 
cases a constant factor of 0.44 was used to go from trans- 
port J 2 x(w)tr to quasiparticle 7 2 x(w)-] A band width of 
W = 1.6 eV was chosen to ensure a zero of Ei(w) around 
~ 400 meV. The results are very similar as to frequency 
and temperature variation to those presented in the top 
frame of Fig. [51 We see, again a resonance like structure 
at ~ 72 meV which decays with increasing temperature 
into a structureless distribution for the real part of the 
quasiparticle self energy Si(w) vs to at 300 K. It is impor- 
tant to note that the resonance peak is seen even above 
T c ~ 91 K in both experiments. [A comparison of the 
superconducting state results for T — 27 K (black solid 
line) and the corresponding normal state results (gray 
solid line) reveals that here the resonance peak is at ~ 62 
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FIG. 5: Top frame: The self energy vs u calculated 

from the spectral densities shown in Fig. U which have been 
found from inversion of ARPES data reported by Zhang et 
all A finite band width of W = 1.2 eV was applied. 
Bottom frame: The self energy Ei(w) vs u calculated 
from electron-boson spectral densities I 2 x(&) obtained from 
optics^ and scaled by a factor of 0.44. A finite band width 
W = 1.6 eV was applied. 



meV. The gap-edge Ao was found to be ~ 19.2 meV and 
is responsible for the shift between normal and super- 
conducting state. For a pure s-wave gap Ao we would 
expect this to be Ao but for d-wave it is less because of 
the distribution in gap values.] 

Finally, we note that in the optics derived case the 
peaks in the lower temperature curves (72 K and 27 K) 
which are in the superconducting state are more pro- 
nounced than they are in the nodal direction ARPES 
data but as we have remarked already, this arises because 
optics is not momentum resolved. We expect that the 
coupling to the optical resonance at ~ 60meV is larger 
in the antinodal direction and that our results are likely 
to be more representative of antinodal direction ARPES 
data. 
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FIG. 6: [Color online] The renormalized dispersion curves 
Ek im meV as a function of the bare energy for various 
temperatures as labeled. The calculations are based in the 
spectral densities / 2 x(w)tr obtained from optics by Hwang et 
aZ.— and scaled down by a factor of 0.44. The band with is 
W — 1.6 eV. Note that in the superconducting state (T = 
72 K and T — 27 K both, nodal and antinodal direction are 
shown. 



In Fig. [5] we compare results for the renormalized en- 
ergy E k vs the bare energy s k using the results of our 
Eliashberg equation solutions based on optics. With 
W = 1.6 eV the crossing is at about 400 meV. For small 
£k the renormalized energy is smaller than its bare value 
which corresponds to an increase in effective mass caused 
by interactions with the medium. This implies band nar- 
rowing. By contrast beyond the zero crossing the renor- 
malized energy is lower than its bare value which corre- 
sponds to band widening. For the infinite band case the 
quasiparticle self energy never crosses zero and becomes 
small only at u) — > oo where bare and interacting disper- 
sion curves meet. There is no concept of band broadening 
or narrowing in this case. We wish to point out another 
interesting feature of these curves. For the two low tem- 
perature curves (solid lines and dashed lines) we are in 
the superconducting state and the full Eliashberg equa- 
tions (|A1[) have been solved with the (i-wave symmetry 
for the superconducting gap built in. In our formulation, 
the self energy Y,(oj + iS) is isotropic, although its value 
does change from its normal state value as the gap opens. 
It is important to understand that for a superconductor, 
the quasiparticle energy is given by 




(8) 



where A^E^) and Z 1 (E] l ) are the real parts of the 
pairing function and of the renormalization function 

ujZ(uj) = Q(oj). Thus, at e k = 0, E k = A(£ k )/Zi(£; k ) 

This is zero in the nodal direction as we can see in the 
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FIG. 7: This illustrates the difference between the self energy 
Ei(u) and the quantity £k — Ek(u) which are different in the 
superconducting state, except in the nodal direction. 



energy in Bi2212 we describe a procedure for extract- 
ing from this information an electron-boson spectral den- 
sity I 2 x(w). This work is based on Eliashberg formalism 
written for a d-wave superconducting gap including finite 
band effects. (The symbol W is used for the band width.) 
The superconducting state is needed because much of the 
data available is for temperatures below T c . Finite bands 
are needed in order to get renormalized and bare bands 
to cross. (They would not in an infinite band formalism.) 
Furthermore, available LDA calculations as well as fits to 
Fermi surfaces also indicate bands of widths of the order 
of some eV. In such cases the final value of the electron- 
boson spectral density I 2 x(^>) from fits to the data is 
significantly affected by W. While we use a maximum 
entropy inversion of the convolution integral (|6|) to get a 
first numerical model for the spectral function, in both 
normal and superconducting state, the Eliashberg equa- 
tions of Appendix [A] are employed with a parameterized 
model for 7 2 x(w) and a least squares fit to the data. 



two curves of Fig. [5] which are labeled 'nodal'. For the 
antinodal direction E^ at £k = is no longer zero. (See 
the two curves in Fig. [5] labeled 'antinodal'.) Finally, we 
note that the quasiparticle lifetime T(Ek) in Eliashberg 
theory is given by 



T(E U 



E k Z 2 (E k ) Ai(£ k )A 2 (£ u ) 



E k Z 2 (E k ) 



(9) 



and this is not simply the imaginary part of the quasi- 
particle self energy in the superconducting state which 
would be E\ L Z 2 {E\ i ) / Z 1 {E] i ) . Because the paring func- 
tion is complex the second term in Eq. @ is non zero and 
this also contributes to the quasiparticle lifetime except 
in the nodal direction for which A^C^k) = 0. While 
Eq. (|A2j) determines the quasiparticle self energy S(w) 
in the normal as well as the superconducting state, the 
equation £i (£?(,/&) = — Ej£ can only be used to deter- 
mine the real part of the quasiparticle self energy in the 
normal state and for the nodal direction in the supercon- 
ducting state. This point is emphasized in Fig. [7] for the 
optimally doped Bi2212 sample BI96A of Hwang et al~ 
with a band width W = 1.6 eV. We compare in this figure 
£i(ct>) and the difference £k — Ek(uj) vs to. In the limit 
lj — > the Si (w) curves go to zero, but £k — Ek(u>) does 
not in the antinodal direction because there is a finite 
gap. This gap was already noted in the data of Kordyuk 
et ali^ The differences are largest at low temperatures. 
In all cases the two curves merge at uj > 100 meV so that 
in this region the difference e^ — E^u;) does not represent 
accurately the self energy in the superconducting state. 



IV. SUMMARY AND CONCLUSIONS 

Motivated by the appearance of new high precision 
ARPES data on the real part of the quasiparticle self 



In making comparison with spectral densities obtained 
from the optical data it is necessary to recognize that op- 
tics involves a momentum average while the ARPES is for 
a single momentum direction (the nodal direction in our 
case). Also, the transport 7 2 x(w)t r can be different from 
its quasiparticle counterpart. In making such a compar- 
ison we noted two main differences. One, the transport 
spectral density is larger in absolute values by a factor 
of about two. Second, the resonant peak around 60 meV 
seen in both spectral densities is more pronounced in op- 
tics. We believe this to reflect the importance of the 
antinodal region not probed in the ARPES data of Zhang 
et al£ Other than these differences there is considerable 
agreement between the two sets of experimental data giv- 
ing some evidence that an Eliashberg approach can be 
used as a phenomenological approach to correlate vari- 
ous data sets. 



As is widely done in this field we interpreted the 
ARPES data directly as the electron spectral function 
assuming the so called 'matrix element effects' described 
by Lindroos et aL— to cause no serious distortion of 
the spectrum. Furthermore, large inhomogeneities are 
seen in the scanning tunneling microscopy (STM) of the 
Bi22I2 compounds. 26 Although this technique probes 
only the surface layer, the inhomogeneities could persist 
in the bulk. Optics is a bulk probe and would average 
over the inhomogeneities so that our analysis would re- 
flect the average spectral density. The good agreement 
with ARPES seen here could be taken as evidence that 
such inhomogeneities, if important, can be treated in an 
average way. An important contribution to the debate 
about mechanism of superconductivity is our finding that 
the coupling to a high energy background seen in optics 
is now confirmed in the ARPES data by Zhang et alZ 
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APPENDIX A: FINITE BAND d-WAVE 
ELIASHBERG EQUATIONS 



for the renormalized pairing potential A(iu) n ,<f>) = 
A( iuj m) cos(20). Here, uj n = nT(2n+l),n = 0, ±1, ±2, ... 
are the electron Matsubara frequencies, T is the temper- 
ature, (...)<£ denotes the average over the polar angle <f> of 
the two dimensional CuO Brillouin zone. Furthermore, 
we have 

oo 

A(m - n) = 2 fdu "f x{y) , 2 (Ale) 

j + {uj m - u n y 





The generalization to d-wave has already been pub- 
lished by Jiang et al^ and has been used to describe var- 
ious aspects of the superconducting state in the cuprates. 
Here we include finite bands as well. Assuming particle- 
hole symmetry for simplicity, the pairs of coupled, non- 
linear, clean limit Eliashberg equations take on the fol- 
lowing form in an imaginary axis notation)^ 



>{iuj n ) = <jj n -\-T A(m— n) 



'u> 2 {iuj m ) + A 2 (iu m , ft) I 
(Ala) 

for the renormalized frequencies Cj(ioj n ) and 



K{iu) n ,4>) — gT ^ cos(20)A(m — n) 



(Alb) 



20(iw m , <j>')A(ioj m , cj)') cos(20') 



Cj 2 (iuj m ) + A 2 (iu; m , <fj) 



and for half-filling 



tan 



W 



2 J Cj 2 (iuj n ) + A 2 (iuj n , , 



(Aid) 



We assumed here, for simplicity, the same form of I 2 x(w) 
holds for the d and the A channel, Eqs. (|Ala|) and (lAlbj) . 
respectively, and the numerical factor g was introduced 
to account for the fact that the projection of the general 
electron-boson spectral density will in general be different 
in the two channels. This factor g can be determined 
from the linearized Eqs. (|Ala[) and (|Alb|) which are valid 
at T = T c whenever T c and I 2 x{uj) are known. In our 
particular case studied here g ~ 1. 

These equations are then to be analytically continued 
using a method formulated by Marsiglio et al£^ and we 
get the following result on the real cj axis: 



~, \ . . T f ,w • x w , ■ m/ 29(iuj m ,(j)')Cj(iuj m ) 

U)(U>) = UJ + 11 y, |A(CJ — lljJm) ~ A(W + VjJ m )\ 



^Cj 2 (iuj m ) + 6 2 {iuj m , (/)') , 



+i J dzI 2 X (z)[n(z)- f(z-io)} 

— OO 

/ 26(u - z)Cj(uj - z + iS) \ 

\ \Jcj 2 ((jj - z + iS) - A 2 (lu - z + iS, 4>') I ; 

for the fully renormalized frequencies Q(u>) and 

x , ,x , , U w ■ ,, u ,. ,J 20(iu m ,cj ) ')A{iLu m ^')cos{2cj ) ') 
A(oj, (j)) = gT cos(20) [A(w - iuj m ) + \{lo + iuj m )\ ( — 

m=o \ ^Cj 2 (iuj m ) + A 2 (iuj m ,<f>') 

oo 

+ig J dz cos(2(j))I 2 x(z) [n(z) — f(z — u)] 

— oo 

/ 29{uj - z)A(uj - z + iS, ft) cos(20') \ 



Jcj 2 (uj - z + iS) + A 2 {uj - z + id, <p') I 

V ch' 



(Ale) 



(Alf) 
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for the fully renormalized pairing potential A(w, </>) in the 
real axis. Here, the function 9(ui,(f>) is defined as 



1(uj, 4>) — tan 1 



and 



AH 



iW 



2JCo 2 {u) - A 2 (w, 



(Alg) 



dv 



I 2 x{v) 
v — oj + id 



(Alh) 



For the normal state only Eq. (jAlaft remains with 
A(iu n , </)) = which is then analytically continued to the 
real axis again using Eq. (|Ale|) with A(u),<j>) set equal to 
zero. The quasiparticle self energy is calculated using the 
relation: 



EH 



LU 



(A2) 



Note that in cases where the square-root is complex, the 
branch with positive imaginary part is to be chosen. 



which is valid in the normal as well as in the supercon- 
ducting state. 
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